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A  new  technique  Is  developed  for  coproxlnatlng  cathode  current-voltage  (I-V) 
and  current-temperature  (I-T)  characteristics  In  a  planar  diode  configura¬ 
tion.  For  a  uniform  work  function  surface,  this  method  allows  both  for 
space-charge  effects  and  for  Schottky  barrier  lowering.  Using  this  method, 
we  derive  an  explicit  first  correction  to  the  Schottky  effect: 
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where  A*  is  the  standard  Schottky  harrier  lowering  (eV),  is  the  net 

enls8lon**current  density,  and  la  the  Child's  Law  current  density.  The 
following  closed  fora  expression  then  approximates  the  net  emission  current 
density: 

"^net  “  “^C* 

•^net  •  •*P  ( V'^c)' 

where  J|^  Is  the  Richardson's  current  density.  This  equation  for  •^net 
connects  the  space-charge  limited  (SCL)  emission  regime  to  the  deep 
temperature-limited  (TL)  regime.  The  SCL-to-TL  transition  point  is  found  to 
be  J{(,  and  the  emission  current,  above.  Is  continuous  with  a  continuous 
slope  at  the  transition  point. 
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I.  INTRODUCTION 


Understanding  cathode  behavior  la  one  of  the  most  crucial  areas  con¬ 
fronting  manufacturers  of  travelling  wave  tube  (TWT)  amplifiers.  The  need  for 
long  lifetime,  high  emission  current,  and  high  reliability  for  potential  space 
applications  means  that  more  stringent  testing  and  evaluation  procedures  are 
needed  to  help  select  suitable  devices.  Here  we  address  one  such  crucial 
area,  that  of  the  expected  cathode  emission  current.  The  net  emission  current 
under  given  specific  conditions  Is  the  single  most  used  criterion  for 
evaluating  cathode  performance.  Until  recently,  full  understanding  of  the 
cathode  I-V  and  I-T  characteristics  for  even  a  uniform  work  function  surface 
was  not  yet  developed.  Expressions  were  known  In  the  limit  of  low  applied 
voltage  or  high  temperature,  and  In  the  limit  of  high  applied  voltage  or  low 
operating  temperature. 

A  full  solution  to  the  uniform  work  function  planar  diode  requires  both 
the  effects  of  space-charge  and  Schottky  barrier-lowering  to  be  Included.  If 
one  Is  to  cover  the  entire  spectrum  from  space-charge  limited  (SCL)  to  temper¬ 
ature-limited  (TL)  emission.  ^Our  work  builds  on  the  historical  worii  of 
Richardson.^  Dushman.^  Child. \  Langmuir.^  and  Schottky.^  and  It  addresses  the 
problem  of  extending  and  comblnxra  the  Langmulr-Chlld  (SCL)  currents  with  the 
Rlchardson-Schottky  (TL)  currents\  It  Is  similar  In  spirit  to  the  work  of 
Nottingham.^  Van  der  Zell.^  CrowelK®  Rlttner^5--****r  thdifttly,  Longo^® 

Scott^^.  and  Hasker.^^ 

We  use  the  same  basic  derivation  method  as  Langmuir^  for  the  space-charge 
effects,  but  we  outline  a  new  technique  for  approximating  the  solution  to  the 
difficult  equations  when  Schottky  barrier  lowering  Is  added.  Here,  an 
approximate  potential  Is  used  to  help  decouple  the  nonlinearity  of  Poisson's 
equation,  which  arises  due  to  space-charge  effects.  We  call  this  method  an 
"effective  potential"  approximation,  since  It  replaces  the  true  potential  with 
an  "effective"  potential  In  one  step  of  the  calculation. 

The  effective  potential  approximation  sacrifices  the  exactness  of  a 
numerical  solution  In  order  to  gain  an  easier  applicability  to  understanding 


emission  data.  In  addition,  an  approximate  analytical  answer  can  be  used  to 
help  make  the  underlying  physical  processes  affecting  the  emission  current 
more  transparent,  and  to  qualitatively  predict  trends  as  a  particular 
parameter  Is  varied. 

In  Sec.  II  of  this  report,  we  outline  some  of  the  physical  assumptions 
made  In  the  present  solution.  In  Sec.  Ill  we  construct  the  approximate 
solution  to  the  nonlinear  Poisson's  equation.  Here  particular  attention  Is 
given  to  how  the  boundary  conditions  should  be  handled.  Indeed,  much  of  the 
simplicity  of  our  final  answer  results  from  carefully  preserving  both  the 
boundary  conditions  and  conserved  parameters. 

In  Sec.  IV,  the  case  of  the  uniform  work  function  surface  Is  considered 
In  detail,  using  this  method.  Me  analyse  how  the  present  technique  reproduces 
the  Child's  Law  current  at  low  applied  fields,  and  how  It  converges  to  the 
Schottky  barrier-lowering  effect  In  very  high  applied  fields. 

The  major  result  of  this  report  Is  then  derived.  A  simple  closed  form 
solution  Is  presented  for  approximating  the  net  emission  current,  which  goes 
smoothly  from  space-charge  limited  (SCL)  to  temperature-limited  (TL)  opera¬ 
tion.  This  new  solution  primarily  corrects  the  Schottky  barrier  lowering 
effect,  to  Include  the  persistence  of  space-charge  effects,  even  at  high 
fields.  This  result  also  shows  that  the  SCL-to-TL  transition  point  Is  the 
Richardson's  current,  and  It  specifies  the  '’roundedness*'  of  the  knee  near  the 
SCL-to-TL  transition  region. 


II.  THE  PHYSICAL  BASIS  OF  CATHODE  I-V  AND  I-T  CURVES 


A.  POISSON'S  EQUATION 

The  effects  of  a  space-charge  cloud  on  cathode  current-voltage  (I-V)  and 
current-temperature  (I-T)  characteristics  are  usually  examined  by  using 
Poisson's  equation: 


V  -  -  p/e^.  (1) 

This  equation  relates  the  local  net  charge  density,  p ,  In  the  space-charge 
cloud,  to  the  electrostatic  potential,  V,  where  Is  the  permittivity  of  free 
space.  The  electrostatic  potential  can  be  converted  to  electrostatic  poten¬ 
tial  energy,  6^,  by  6^  >  q^  V,  inhere  q^^  is  the  charge  of  Che  particle  under 
the  Influence  of  V.  Also,  p(r)  for  the  surrounding  medluni  can  be  written  as: 

p(r)  •  q^  n(r),  (2) 

where  n(r)  Is  the  positive  number  density  of  particles  In  the  medium,  and  qj 
Is  .their  Individual  charge.  For  electrons  interacting  with  an  electron  space- 
charge  cloud,  Poisson's  equation  becomes: 

V^0(r)  ■  -e^  n(r)/e  ,  (3) 

o 

—  IQ 

where  e  ■  +1.602  x  10  coulombs. 

When  Image  force  effects  are  Included,  n(r)  will  contain  both  a  bulk 
term  n(r),  describing  Che  average  properties  of  the  space-charge  cloud,  and  a 
singular  term  6(r),  which  Is  a  Dirac  delta-function  term  describing  the 
electrode  reaction  to  an  Individual  discrete  electron. 

An  equivalent  formulation  Is  to  first  solve  Poisson's  equation  In  the 
zero  current  case  for  0  where: 


and  where  p  (r)  contains  all  the  singular  aspects  of  the  charge  dlstrlbu- 
tlon.  The  monopole  contribution  to  image  forces  on  an 

individual  electron,  and  to  charges  on  the  plates  from  placing  an  external 
applied  voltage  between  the  electrodes.  The  dipole  contribution  to 
due  to  the  surface  work  function  (^),  or  work  function  distribution  on  the 
electrodes. 

When  there  is  a  net  current  flow  between  the  plates,  one  can  then  solve: 

[e(r)  -  0^^^  (r)l  -  -e^  n(r)/e^,  (5) 

where  e(r)  is  the  total  electron  potential  energy  at  any  point  r  between  the 
electrodes,  and  n(r)  is  the  local  bulk  space-charge  density  at  that  point, 
which  is  non-singular.  Using  the  0._,.(r)  function  assures  convergence  to  the 
proper  solution  when  the  bulk  space-charge  density  n(r)  vanishes.  The  n(r) 
bulk  charge  density  would  vanish,  for  example,  when  the  cathode  temperature  is 
lowered  under  constant  applied  field  conditions,  as  in  a  dip  test. 

We  next  simplify  the  problem  to  one  rectilinear  dimension,  x,  which  is 
the  planar  diode  configuration,  resulting  in: 

j2  2 

— j  [0(x)  -  n(x)/e^.  (6) 

dx^ 

Often,  there  is  only  one  direction  of  interest,  defined  by  the  average  net 
current  flow  direction,  but  in  many  applications  the  electrodes  may  contain 
curved  surfaces.  However,  since  much  of  the  net  current  determination 
Involves  what  happens  within  microns  of  the  cathode  surface,  using  an  "effec¬ 
tive  distance"  parameter  in  Eq.  (6)  may  help  the  planar  diode  solution  to  be 
applied  to  these  curved  geometries. 

B.  APPROXIMATIONS  USED  IN  DERIVING  THE  SOLUTION 

What  makes  the  solution  to  Eq.  (6)  difficult  is  that  the  local  charge 
density  between  the  plates,  n(x),  is  determined  by  the  amount  of  current 
flowing,  which  in  turn  is  determined  by  various  aspects  of  the  0 (x)  curve 


Itself.  This  effect  change's  n(x)  Into  n(6(x)),  naklng  Eq.  (6)  Inherently 
nonlinear.  In  the  effective  potential  approximation  technique,  this 
nonlinearity  is  circumvented  by  replacing  the  6(x)  function  in  the  n(9(x)) 
term  with  an  effective  potential  function  which  retains  the  most  important 
characteristics  of  the  true  0(x)  function.  This  technique  is  much  like 
replacing  a  known  non-elementary  function  by  a  polynomlnal  fit  to  that 
function. 

A  second  difficulty  in  Eq.  (6)  is  that  the  d(x)  function  necessarily  has 
a  maximum  in  it,  at  some  point  x^  between  the  electrode  plates.  The  maximum 
point  in  0(x).  denoted 

max  t0<x)l  5  e(x^)  5  E^,  (7) 

then  represents  a  potential  energy  barrier  which  electrons  must  overcome  in 
order  to  be  registered  as  collected  current.  Electrons  emitted  tflth  an  energy 
lower  than  this  barrier  height  E^  will  be  reflected  back  into  the  cathode. 

The  reflected  electrons  do  not  contribute  to  the  collected  current,  but  they 
can  alter  the  local  n(x)  net  electron  density,  forming  a  "trapped  cloud"  in 
the  near  cathode  region. 

In  much  of  this  work,  we  neglect  the  contribution  to  n(x)  due  to  elec¬ 
trons  in  the  trapped  cloud.  It  is  equivalent  to  assuming  that  only  electrons 
with  enough  energy  to  surmount  the  barrier  are  emitted.  He  shall  see  later 
that  at  low  voltages  the  present  calculation  for  x^  and  E^  can  give  x^  values 
which  can  be  a  significant  fraction  of  the  total  anode-to-cathode  distance 
d.  Such  a  wide  extent  for  the  barrier  position  means  that  if  the  calculation 
Included  the  effect  of  a  few  electrons  forming  a  trapped  cloud  with  energies 
near  E^,  they  would  not  change  the  barrier  height  by  much.  However,  by  inter¬ 
acting  among  themselves,  they  could  alter  the  barrier  position  signifi¬ 
cantly.  Thus,  by  neglecting  these  electrons  forming  a  trapped  cloud,  we 
expect  to  overestimate  the  barrier  position  x^  by  a  fair  margin,  especially  in 
the  space-charge  limited  region,  but  the  net  barrier  height  E^  may  not  be 
affected  as  much. 


III.  CONSTRUCTION  OF  THE  APPROXIMATE  SOLUTION 


A.  THE  ELECTRON  DENSITY  IN  THE  SPACE-CHARGE  CLOUD 

Figure  1  shows  the  potential  energy  diagram  for  electrons  between  the  two 
electrodes.  The  zero  of  potential  energy,  E-O,  Is  referenced  to  the  vacuum 
level,  l.e.,  the  potential  energy  at  large  distances,  extrapolated  to  zero  net 
field  between  the  electrodes.  The  x«0  point  Is  the  cathode  emitter  plate,  and 
x-d  Is  the  anode  collector  plate.  The  net  barrier  height  and  the  barrier 
position  are  labelled  E^.  and  x^,  respectively.  We  next  develop  the  approxima¬ 
tion  techniques  to  estimate  both  E^  and  x^,  given  a  set  of  conditions,  such  as 
emitter  work  function  temperature  T,  applied  voltage  Va»  and  anode-cathode 
distance  d. 

Poisson's  equation  [Eq.  (6)],  must  be  solved  In  the  presence  of  two 
conserved  quantities: 

a)  Conserved  total  electron  energy: 

j  m  v(x)^  +  6(x)  -  E,  (8a) 

v(x,E)  -  {2(E-e(x)J/m}^/2,  E>  0  (x).  (8b) 

b)  Conserved  forward  current.  Independent  of  position: 

dJ  -  (A^/k)  T  exp  I-(E-t+^)/kT]  dE,  (9a) 

+  o  e 


^  - 


Awmek 


s  119.58  amps/(cm^  K^), 


(9b) 


where  A^  Is  Richardson's  constant  and  k  Is  Boltzmann's  constant.  The  net 
current  Is  derived  from  Eq.  (9a)  by  summing  all  the  current  contributions  from 
E  -  (E^,,  -). 
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ZERO  CURRENT 
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Schematic  of  the  inter-electrode  potential  energy  vs.  distance,  zero 
current,  and  finite  current  cases.  The  emltter^and ^collector  work 
functions  («  and  «  )  and  their  Fermi  levels  (E«.  sp  are  shown, 
along  with  tfie  applied  voltage  V^.  The  space-cnarge  associated  with 
a  finite  emitter- to-collector  current  causes  an  electron  barrier  to 
develop,  denoted  E^.  which  can  be  several  microns  away  from  the 
cathode,  at  a  position  labelled  x^.  Only  electrons  emitted  with 
energies  greater  than  E^  are  collected  as  the  emitted  current 


density. 
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In  Poisson's  equation,  the  local  number  density  of  electrons  n(x)  can  be 
written  as: 

n(x)  -  /  n(x,E)  dE.  (10) 

(all  E) 

The  n(x,E)  term  Includes  all  electrons  with  an  energy  greater  than  E^,  which 
become  part  of  the  collected  current,  as  well  as  those  electrons  with  energies 
less  than  E^,  which  could  form  a  trapped  cloud  near  the  emitter.  Neglecting 
the  contribution  to  n(x)  due  to  electrons  with  energies  E  <  E^  constitutes  the 
"empty  cloud"  approximation  used  here.  Equation  (10)  then  becomes: 


n(x)  •  /  n  (x,E)  dE, 

E  -  E 


(11a) 


n  (x,E)  e  v(x,E)  dE 
n  (x.E)  dE - .  - 


e  V  (x,E)  • 


(11b) 


Direct  substitution  of  Eqs.  (8)  and  (9)  Into  Eqs.  (lla)>(llb)  yields: 


2A  21/2 

n<x)  -  (i  (-2^ - )  (^)  f(y/)  exp  HE^  ♦  ♦^)kT] ,  (12a) 

f(y  ^)  s  (4/w)^^^  exp  (+  y  /  dy  exp  (-y^)  (12b) 

y-Xg 

2 

-  exp  (+  Vf,  )  erfc  (y^), 

y^2(x)  5  [E^-8(x)]/kT.  (12c) 

Defining: 

1/2 

"o  “  [rcrl  , 


(13) 


it  then  corresponds  to  the  density  of  electrons  at  the  barrier  maximum,  with  c 
being  the  speed  of  light.  Equation  (13)  also  shows  that  this  electron  density 
Is  proportional  to  the  net  collected  current.  Poisson's  equation,  using  Eq. 

(I la),  then  becomes: 

-e^n  y 

(e(x)  -  e^^^(x)i  -  fty^^x)].  (14) 

dx  o 

The  approximate  solution  to  this  equation  is  developed  next. 

B.  SOLUTION  TO  POISSON'S  EQUATI(»1  AND  BOUNDABY  CONDITIONS 
The  solution  to  Eq.  (14)  can  be  formally  be  written  as: 

»(.)  -•e.„<x)  -  V  ♦  'o  - 

e^n  XX.  , 

— —  /  dx.  /  *  dx,  f[y  ^(x,)],  (15) 

*o  Xj-0  ^  Xj-O  2 

where  and  B^  are  constants.  The  double  Integral  can  be  transformed  Into  a 
single  Integral 

9(x)  -  0_(x)  -  K  X  +  B^  - 
ext  o  o 

2 

e  n  X 

— —  /  (x-x)  dz  f(y  ^(s)).  (16) 

o  x«0 

The  KqX  and  B^  terms  In  Eq.  (16)  give  the  electrode  reaction  to  the  presence 
of  the  space-charge  cloud  between  the  plates. 

The  boundary  conditions  at  the  electrode  plates  are  next  applied.  To 
preserve  0(x“O)  •  must  set  Bq*0.  To  preserve  0(x»d) 

■  requires  that: 

«^n  d  - 

/  (d-z)  dr  fly  ^x)l,  (17) 

o  r“0 

determining  the  value  for  K^. 
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Next,  define  two  helping  functions: 


g(x)  -  /  fty_  (z)]  dzt 
z-0 

X  - 

h(x)  -  /  f[y  (*)]  t  dz, 
z-0 


(18a) 

(18b) 


and  one  can  derive  the  following  expressions  for  6(x)  and  d6(x)/dx: 

4e^n 

e(x)  -  0g^t(x)  “  ~nr  ^  t8(d)-g(x)]  +  h(x)  d  -  h(d)  * 


de(x) 

dx 


d0  ,(x) 
ext 

dx 


+e  n 

{d  [g(d)  -  g(x)l  -  h(d)}  . 


(19a) 


(19b) 


From  these  two  equations,  two  other  combinations  can  be  constructed.  The 
first  combination  is: 


Atk  f  *  ® 

[e(«)  -  -  X  ^  (h(,)) .  (20) 

o 

and  it  Involves  the  nonlinear  part  of  Poisson's  equation  only  In  the  region 
between  x-0  and  the  present  value  of  x,  through  the  h(x)  function.  The  second 
combination  Is: 

(6<x)  -  ♦  (J-x)  (f  -  fgl)  . 

+e^n 

(g(d)-g(x)l  -  [h(d)-h(x)l}  ,  (21) 

and  this  Involves  the  nonlinear  part,  through  the  g(x)  and  h(x)  functions, 
only  in  the  region  between  the  present  value  of  x  and  the  collector  at  x-d. 

The  general  shape  of  6(x),  as  shown  In  Fig.  1,  has  a  single  maximum 
within  the  range  x-(0,d),  occurring  at  x  -  x^,  with  max[6(x)]  =  E,..  Since 
Eqs.  (20)  and  (21)  are  valid  for  all  x.  It  is  more  convenient  to  use  Eq.  (20) 


for  describing  the  region  of  x-(0,x^),  and  to  use  Eq.  (21)  for  describing  the 
region  x  ■  (x^.d).  By  deriving  equations  which  separate  the  region  x>(0,Xg) 
froo  the  region  x-(Xg,d),  one  can  easily  apply  the  condition  which  requires 
x*Xg  to  correspond  to  the  maximum  point  In  6(x). 

The  separation  Into  x<x^  and  x>x^  regions  also  allows  the  solution  to  be 
expressed  In  terms  of  scaled  distance  variables.  It  makes  the  evaluation  of 
expressions  In  the  critical  x  x^  region  easier.  We  use  the  following  scaled 
variables : 


Y  -  x/d,  Yq  -  x^/d, 

(22a.b) 

U  -  (x/x^)  -  1, 

(22c) 

W  -  I(x-x^)/(d-x^)] , 

(22d) 

so  that  0  •  (-1,  0)  when  x  •  (0,  x^),  and  W  «  (0,  +1)  when  x  - 

(Xg.d).  The 

terms  in  Eqs.  (20)  and  (21)  Involving  the  nonlinear  functions, 
can  then  be  expressed  in  terms  of  the  following  functions: 

g(x)  and  h(x). 

,  U 

Gi(u  -  i-  -  l)  -  /  f[ul  du. 

c  u— 1 

(23a) 

,  U 

Hi(u  -  1)  -  /  f[ul  u  du. 

C  U"-l 

(23b) 

■  3^)  ■  i,’ 

c  w»W 

(23c) 

x-x  1 

H2(W  -  ■  /  f  [wj  w  dw. 

c  w-W 

(23d) 

d 


By  using  these  scaled-variable  functions.  Eqs.  (20)  and  (21)  become: 


[«(-)  -  (.)1  -  t  d  (H  -  ^) 


e^n 

(— r - )  [Gj(U)  +  Hj(U)]. 


(24a) 


de 


^  d^ 

-  (^-^2 - )  (1  - 

o 


(24b) 


The  last  boundary  condition  can  now  be  easily  applied, 
at  X  Xgt  both  the  following  conditions  oust  be  satisfied: 


It  requires  that 


(25a) 


9(x)  ■  (at  X  «  x^). 
c  c 


(25b) 


Applying  these  conditions  to  Eqs.  (24a)  and  (24b)  then  results  In  the  follow¬ 
ing  two  algebraic  equations,  for  the  two  unknowns,  E^,  and  x^: 


[e  -e  ^(x  )]-  Y  (-d  T-^ 

c  ext  O  '*  o  '  dx  'x 


aV  d^  , 

(-j2-)  Tr„  [0,(0)  +  H,(0)] 
o 


(26a) 


de  .  I 

♦  (•  -  i-o)  H  3^  I. ) 


-  (— T— )  (1  [G2(0)  -  H2(0)]. 

o 


(26b) 


This  method  has  manipulated  the  nonlinear  Poisson's  equation  to  emphasize 
the  two  parameters  of  interest:  the  net  barrier  height  E^,  and  the  barrier 
position  x^.  It  also  shows  what  part  of  the  full  solution  to  the  differential 
equation  couples  into  the  E^  and  determination.  Equations  (26a)  and  (26b) 
show  that  the  nonlinearities  involve  only  the  [6^(0)  Hj^(O)]  and  [62(0)  - 

H2(0)]  terms,  both  of  which  are  independent  of  y ^  -  x^/d.  What  we  have  done 
is  to  preserve  the  boundary  conditions  and  symmetries  Inherent  in  the  original 
problem,  while  seeking  out  E^  and  x^. 

C,  THE  EFFECTIVE  POTENTIAL  APPROXIMATION 

Equations  (26a)  and  (26b)  are  now  in  a  suitable  form  for  using  an  effec¬ 
tive  potential  approximation.  Since  the  nonlinear  terms  (Gp  G2,  Hj>H2) 
contain  only  E^  and  Integrals  over  whole  portions  of  the  6 (x)  function,  the 
evaluation  of  Gj(0),  H^(0),  G2(0),  and  H2(0)  in  Eqs.  (26a)  and  (26b)  is  less 
sensitive  to  the  exact  details  of  the  6(x)  function  than  if  these  terms  con¬ 
tained  terms  involving  6(x)  directly  or  derivatives  of  6(x).  The  potential 
energy  curve  d(x)  enters  into  the  nonlinear  terms  through  the  definition  of 

yc^(*)  -  “  e(z)]/kT,  (27) 

where  z  is  the  Integrand  variable  in  each  non-linear  function. 

To  evaluate  the  nonlinear  terms  in  Eqs.  (23a)-(23d),  we  approximate  the 
true  6(z)  function  in  these  integrals  by  a  polynomial  in  z,  or  by  a  power-law 
function  in  z.  This  replacement  of  the  true  6(z),  which  is  a  difficult  tran¬ 
scendental  function,  by  an  approximating  elementary  function,  is  equivalent  to 
substituting  an  effective  potential  in  for  6(x).  The  choices  available  for 
the  effective  potential  function  can  be  narrowed  by  requiring  the  function  to 
obey  known  critical  properties  of  the  true  6(z)  function. 

For  the  region  x  -  (0,  x^)  three  boundary  conditions  on  6 (x)  must  be 
satisfied  by  the  effective  potential: 
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(28a) 


§ 


0(x-O)  -  - 


6(x-x  )  -  E  . 
c  c 


(28b) 


de(x) 


lx  'X 


(28c) 


The  simplest  polynomlnal  function  In  x,  which  obeys  these  three  conditions,  is 
the  parabola: 


e(x)  -  -  +  E^)  [1  -  (x/x^)]^  +  E^. 


The  parabola  is  also  special,  because  it  corresponds  to  an  approximately 
uniform  net  charge  density  near  the  cathode,  making  the  near-cathode  region 
much  like  a  “virtual  cathode"  from  which  emission  current  can  be  drawn.  This 


parabolic  approximation  then  gives: 


y^iz)  -  (1  -  <*/x^)J^  (♦^  +  E^)/kT. 


(30a) 


y„^(u]  -  u2  (*  +  E J/kT 


(30b) 


Equation  (30b)  can  then  be  substituted  into  the  and  integrals  [Eqs. 
(23a)-(23b)]. 


For  the  region  x  >  (x^,d),  the  following  boundary  conditions  need  to  be 
satisfied  by  the  effective  potential: 


6(x-x  )  -  E  , 
c  c 


(31a) 


de(x) 

dx 


(31b) 


0(x-d)  = 


-®E--  (^^^♦e-*c>* 


(31c) 


where  Is  the  applied  external  voltage,  and  is  the  collector  or  anode 
work  function.  Here,  we  could  also  use  a  polynomial  to  approximate  the  9 (x) 
function  for  the  x  •  (Xj.,d)  region.  However,  an  additional  piece  of  Informa¬ 
tion  Is  available  from  Child's  Law,  which  says  that  for  large  x,  away  from  the 
emitter: 

e(x)  -  (const.)  x^^^.  (32) 

We  thus  choose  the  following,  as  the  simplest  effective  potential  for  the  x  ■ 
(x^,d)  region,  and  one  which  satisfies  Eqs.  (31a, b,c)  and  (32): 


9(x)  -  -  (0_  +  E  )  l(x-x)/ 

C  El  C  C  C 


This  choice  then  gives: 


y^2(2)  .  [(z-x^)/(d-x^)]^/^  (6g  +  E^)/kT, 


(34a) 


(Sg  +  E^)/kT. 


(34b) 


Equation  (34b)  can  then  be  substituted  Into  the  G2  and  H2  Integrals  [Eqs. 
(23c)-(23d)]. 

Grouping  the  nonlinear  terms  together  In  the  main  equations  [Eqs.  (26a)- 
(26b)],  and  defining  new  functions,  f|^  and  f^: 


5  [Gj(0)  +  Hj(0)], 
f  =  [G2(0)  -  H2(0)], 


(35a) 


(35b) 


gives  the  following  for  the  function  f|^: 


fk[A  -  (♦e  +  E-)/kT]  -  /  (1+u)  du  f(y^  -  Au^), 


1 


r 


and  it  gives  for  fy; 


fy[A  -  (Sg  +  E^)/kT]  -  (1-w)  dw  f(y^  -  Aw^'^). 


Finally,  the  function  f(y-^)  In  the  Integrand  [see  Eq.  (12b)]  can  also  be 


approximated  by  algebraic  functions. 


Abramowltz  and  Stegun^^  place  bounds  on  the  f(y_^)  function: 


2  TTi  ^  *'^c  '  ^  2  TTT  ~  ^o'^c  '* 

Yc  +  (Fc  +  2)^'“^  +  (Yc  + 


The  upper  bound  In  Eq.  (38)  at  y^  =  0  corresponds  to  the  point  0  (x)  ■>  E^,.  We 
use  this  upper  bound  to  the  f (y^^)  function,  denoted  ^^(yj.^) >  here  because  It 
overestimates  both  the  density  of  charges  and  the  magnitude  of  space-charge 
effects.  It  compensates  In  part  for  not  explicitly  considering  electrons 
forming  a  trapped  cloud  near  the  cathode  [see  Eqs.  (lO)-(ll)],  which 
underestimated  the  space-charge  cloud  effects. 


Using  the  fQ(y^^)-f unction  then  gives: 


(39a) 


f  [A  .  (6  +  E  )/kTl  .  /  (1-.)  d» - 4/S  4/3 - iTI- 

^  w-  0  (A  +  [AW^'-*  +  (4/ff  )]^'^ 


(39b) 


as  explicit  Integrals  for  approximating  the  nonlinear  terms  In  Poisson's 
equation.  These  Integrals  are  worked  out  In  the  appendix. 


igggggSgisas  ■ 


D.  DETERMINING  THE  BARRIER  HEIGHT  AND  BARRIER  POSITION 

In  Eqs.  (39a)  and  (39b)  both  fy  and  are  functions  of  the  net  barrier 
height  Ej.,  but  not  of  the  barrier  position  x^.  The  fy  and  fj^  functions 
contain  all  the  effects  of  the  nonlinearity  of  the  original  Poisson's  equa¬ 
tion.  In  addition,  the  fy  and  fj^  functions  naturally  separate  the  solution 
Into  a  cathode  work  function  term,  )>  and  an  applied  field  term, 

fy(E^,  9g).  If  one  were  to  use  a  different  effective  potential  function  than 
those  chosen  here  [see  Eqs.  (29)  and  (33)],  the  fy  and  fj^  functions  would 
still  be  of  the  general  form: 


fy(A)  «  (const. )/A^^^, 


(40a) 


tj^(A)  •  (const. )/A^^^. 


(40b) 


Additional  properties  of  Eqs.  (26a)  and  (26b)  can  be  used  to  help  make 
the  determination  of  the  roots  E^  and  x^.  easier.  Using  the  following 
definitions: 


«o  -  ®c  -  ®ext(*c>' 


(41a) 


„  /  X  .  ext 

-sr- 


(41b) 


e2  d2 

m  ■  - n  , 

o  e  o' 
o 


(41c) 


then  allows  Eqs.  (26a)  and  (26b)  to  be  rewritten  as: 


(42a) 


a  -  Y  B  •  m  Y  i (E  ) . 

O  O  O  O  *0  iC  c 


(42b) 


24 


The  term  is  a  strong  function  of  and  it  is  Independent  of  position: 


m  -  m_  exp  (-  E  /kTj ,  (43a) 

OK  C 

-R  =  (^)  (r-ET)  (¥-)  N.'“)  • 

o 

The  Og  term  is  an  implicit  function  of  both  Xg  and  Eg,  and  it  represents  the 
net  barrier  height  Increase  caused  by  the  presence  of  the  space-charge 
cloud.  For  accelerating  field  conditions  (6^  ^  positive,  resulting 

in  the  following  inequality: 

>  ^  .  (44) 

This  inequality  can  be  used  to  help  narrow  the  range  of  possible  x-values  to 
regions  much  smaller  than  x«*(0,d),  when  seeking  the  roots  of  Eqs.  (42a)  and 
(42b).  Equation  (44)  also  shows  how  the  presence  of  the  space-charge  cloud 
has  naturally  placed  an  upper  bound  on  how  far  Xg  can  tiander  away  from  the 
cathode. 

Solving  for  m^  and  in  Eqs.  (428)  and  (42b),  then  gives: 


(45a) 

(45b) 


Equation  (45a)  is  a  transcendental  equation  which  can  determine  Xg.  It  has  a 
weak  dependence  on  Eg.  due  to  the  nonlinear  space-charge  functions,  fy  and 

For  a  given  value  of  Xg,  Eq.  (45b)  then  gives  the  barrier  height  Eg.  The 
true  roots  of  Eqs.  (45a)  and  (45b),  Eg*  and  Xg*,  can  now  be  easily  determined 
using  the  following  procedure. 
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mi 

Begin  with  and  determine  the  Intrinsic  barrier  maximum.  In  the 

absence  of  any  space-charge  effects: 

max  (e^^^.(x)]  =  (A6) 

Using  as  an  Initial  value,  fy(Eg^)  and  are  evaluated  using  Eqs. 

(39a)-(39b)«  Equation  (45a)  then  determines  the  x^  position  based  on  the 
properties  of  the  original  function  and  Its  first  derivative  (through 

the  term).  This  x^  value  can  then  be  substituted  Into  Eq.  (45b),  along 
with  the  and  £)((Ec^)  values,  to  determine  a  modified  barrier  height 

value,  Eg^.  The  E^^  can  now  be  used  in  place  of  E^^,  and  the  process 
repeated.  After  a  few  Iterations,  this  procedure  converges  to  the  values  E^* 
and  which  are  the  real  roots  of  Eqs.  (45a)  and  (45b). 

E.  RELATIONSHIP  OF  BARRIER  HEIGHT  TO  EMISSION  CURRENT 

Equations  (39a)-(39b)  and  (45a)-(45b)  represent  one  major  result  of  this 
report.  The  values  E^*  and  x^*  for  the  barrier  height  and  position  are 
determined  using  the  particular  effective  potential  approximation  which  Is 
contained  In  the  f^  and  f|^  terms.  To  change  the  effective  potential  only 
alters  the  fy  and  f|^  Integrals,  the  rest  of  the  procedure  remaining  the 
same.  In  this  manner,  the  original  nonlinear  Poisson's  equation  Is  reduced  to 
root-finding.  The  collected  current  density,  which  Is  determined  by  the  total 
amount  of  electrons  which  can  surmount  the  barrier.  Is  then  given  by: 

^net  ” 

where  A^^  Is  Richardson's  constant,  and  Is  the  emitter  work  function.  From 
Eg*  and  Eq.  (47),  approximate  current-voltage  (I-V)  and  current-temperature 
(I-T)  characteristics  can  be  constructed. 


IV.  THE  UNIFORM  WORK  FUNCTION  CASE 


A.  THE  SPACE-CHARGE  AND  TEMPERATURE-LIMITED  CURRENTS 

In  this  section,  the  effective  potential  method  for  generating  approx¬ 
imate  I-V  and  I-T  curves  Is  applied  to  the  uniform  work  function  emitter,  with 
both  Schottky  barrier  lowering  and  space-charge  effects  Included.  We  first 
define  a  Richardson's  Current  by: 

-  A^  t2  exp  [-♦^/kT] .  (48) 

and  we  also  define  a  Child's  Law  Current  Jq  by: 


^  9—  1/9  9/9  ^  CB»/e^ 


3/2 


^  9  d2 

where  K  Is  2.33S2S  x  10"^  amps/volt^^^,  and  where  6^  Is: 


®E  -  •  ^  ^  ♦e  -  ♦c' 


(49) 


(50) 


Is  the  applied  external  potential,  and  4^  and  4^  are  the  emitter  and 
collector  work  functions.  The  and  currents  help  organise  the  parameters 
In  the  effective  potential  approximation  to  more  clearly  exhibit  the  transi¬ 
tion  from  space-charge  limited  (SCL)  to  temperature-limited  (TL)  behaviors. 

The  effective  potential  method  begins  with  specifying  the  potential 
energy  function  for  an  Individual  electron  In  the  limit  of  no  net  current, 
6^^^(x).  For  a  uniform  applied  field  with  a  near-cathode  Image-force 
potential,  c®"  I*®  written  as: 

0  (x)  -  -  (x/d)  0  -  e^/[l6x  e  (x4x  )],  (51) 

ext  E  o  a 
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where  Is  an  "added  distance"  which  prevents  the  Image  force  potential  from 
being  singular  at  the  cathode  (x"0).  The  boundary  condition  at  the  cathode 
sets  the  value  of  x^  via: 

e^/Ilb  x^l  -  (52) 

For  a  work  function  of  -  2  eV,  x^*  1.8  A.  An  additional  image  force  for  ' 
the  near-collector  region  can  also  be  included  in  Eq.  (51)  for  increased 
accuracy  at  near-retarding  field  voltages. 

The  effective  potential  method  presents  two  algebraic  equations.  Eqs. 
(42a)-(42b)  or  Eqs.  (45a)-(45b),  whose  roots,  E^  and  x^,  represent  the 
approximate  barrier  height  and  barrier  position,  in  the  presence  of  space- 
charge  effects.  The  critical  parameters  from  the  external  applied 
potential,  .**■«  1'^®  value  at  x^  and  its  derivative  at  x^..  For  the 

potential  of  Eq.  (51),  the  derivative  term  is: 


d  6 


ext 

3x" 


-®E- 


e^  d 


16  we^  (x  +  X  ) 
o  c  a 


I  • 


(53) 


The  nonlinear  coupling,  i.e.,  how  the  space-charge  affects  the  current 
between  the  cathode  and  the  anode,  was  expressed  in  terms  of  two  additional 
functions,  labelled  fy  and  fj^,  which  both  depend  on  the  barrier  height  E^,  but 
not  the  barrier  distance. 

The  function  fj^  is  only  a  function  of  the  combined  parameter  A  -  (E^ 
4^)/kT,  and  as  we  show  later,  the  essential  transition  from  SCL  to  TL  opera¬ 
tion  does  not  critically  depend  on  fj^.  The  function  fj^  turns  out  to  be  more 
critical  in  the  very  low  field  regime,  where  there  is  another  transition  from 
SCL  operation  to  what  we  will  call  almost-retardlng-field  (ARF)  operation. 

This  ARF  regime  is  characterized  by  significant  deviations  from  a  Child's  Law 
type  behavior. 


The  function  fy  Is  only  a  function  of  the  combined  parameter  A  * 

(Eg-t-  9g)/kT.  For  applied  voltages  greater  than  a  few  volts  (l.e.,  outside 
ARF,  but  allowing  either  SCL  or  TL  regimes),  one  will  automatically  be  in  the 
large  A  limit  of  fy.  As  shotm  in  the  appendix  [Eq.  (A-14)],  for  large  A: 


3  C 

fy  (A)  -  |g; - +  [Order  (r"^)]  ,  (54) 

2  R 

where  R^  =  (xA/4). 

In  Eq.  (54),  the  parameter  is  an  integration  constant  which  depends  on 
the  details  of  the  entire  inter-electrode  potential.  The  effective  potential 
method,  by  using  an  approximate  function  for  the  inter-electrode  potential, 
can  result  in  an  error  in  the  determination  of  C^.  As  such,  Eq.  (54)  may  only 
be  accurate  to  the  first  term.  Itowever,  we  show  next  that  the  leading  term  in 


,  »  (e.  +  H.) 

•  (55> 

is  sufficient  to  generate  a  SCL  to  TL  transition. 

The  transition  from  SCL  to  TL  behavior  can  be  shown  most  easily  using 
Eqs.  (42a)-(42b).  Both  those  equations  are  expressed  in  terms  of  a  scaled 
distance,  where  x^  is  the  distance  to  the  barrier  maximum,  and  d  is 

the  anode-cathode  separation.  For  conditions  that  are  not  ARF,  Is  expected 
to  be  small.  Neglecting  terms  of  order  in  Eq.  (42a)  gives: 


m  (1  -  2  Y  )  f_.. 

(56a) 

6 

o 

x 

m  exp  (-E  /kT). 

K  C 

(56b) 

The  parameter  m^^  [see  Eq.  (43b)]  can  be  expressed  in  terms  of  the  currents  Jj^ 
and  as  follows: 
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(57) 
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Substituting  Eqs.  (53),  (55),  (56b),  and  (57)  Into  Eq«  (56a),  and  neglecting 
X3  as  small  compared  to  x^,  then  gives: 


«E- 


e^  d 


16  ir  €  X  ^ 
o  c 


J»  0-  1/2  2  x^ 

®E  jT  (e,  ^  E  )  "”d“^  [-E^/kT]. 

C  E  c 

The  near  unity  terms  In  Eq.  (58)  can  be  further  approximated  by  unity  for 
examining  the  SCL  to  TL  transition,  giving: 


exp  (-E^/kT)  s 


(58) 


(59) 
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16  X  e  X 
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Equation  (59),  by  itself,  embodies  a  SCL-to-TL  transition.  When  the 
barrier  position  wanders  on  the  order  of  l-nim  away  from  the  cathode,  for  d  ■ 

1-cm,  the  x^  correction  term  in  Eq.  (59)  is  small  (•  3.6  eV).  Under  these 

conditions,  Eq.  (59)  shows  that  "^net  ■  which  Is  SCL  operation.  For  higher 

applied  fields,  the  barrier  maximum  moves  closer  to  the  cathode,  and  the  x^ 
correction  term  In  Eq.  (59)  then  becomes  a  significant  fraction  of  6^.  The 
net  resulting  current  Is  then  less  than  Child's  Law  would  have  predicted  at 
those  applied  voltages,  and  one  Is  In  TL  operation. 

Examining  the  other  equation  Involving  x^  and  E^  [Eq.  (42b)],  it  has  a 
term  In  It  which  Is  solely  a  function  of  the  external  applied  poten¬ 
tial,  9g^t(*)»  Eor  the  function  here  [Eq.  (51)],  the  terms  in  Eq, 

(42b)  can  be  combined  as  follows: 


«.,t  'V  »„■-  2*  '.«c). 


(60) 
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which  gives  for  Eq.  (42b): 


-  -  t2e^/(16  ir  x^)l  +  (x^/d)^  exp  (-E^/kT).  (61) 

The  contribution  to  the  net  barrier  height  E^  in  Eq.  (61)  from  Eq.  (60)  is 
necessarily  negative,  since  it  corresponds  to  Schottky  barrier  lowering.  The 
contribution  to  E^  involving  fj^  is  necessarily  positive,  and  corresponds  to 
raising  of  the  electron  barrier,  due  to  space-charge  effects.  Similar  to  how 
Eq.  (59)  changes  %flth  x^  value,  one  sees  that  large  x^  values  give  a  large 
space-charge  term  and  a  small  Schottky  barrier  term.  This  makes  E^  tend  to  be 
positive  giving  rise  to  SCL  emission.  Conversely,  small  x^  values  are  asso¬ 
ciated  with  a  small  space-charge  term  and  a  large  Schottky  barrier  term,  cor¬ 
responding  to  TL  eoflsslon. 

For  SCL  emission,  E^  can  be  positive.  Equation  (61)  then  shows  that  the 
fj^-term  has  a  central  role  in  allowing  E^  to  be  positive.  In  contrast,  Eq. 
(59),  and  its  "parent"  equation,  Eq.  (42a),  have  only  a  weak  or  no  dependence 
on  f^^.  Thus,  when  Eq.  (59)  reduces  to  Jn^t  “  ‘^C  regime,  it  means 

that  Eq.  (59)  is  practically  determining  E^  by  Itself,  since  *^net  contains 
E^.  With  E^  almost  exactly  specified  in  SCL  emission  by  Eq.  (60),  the  role  of 
Eq.  (61)  and  thus  the  role  of  the  fj^-function  in  SCL-emlssion  is  then  changed 
into  primarily  determining  barrier  position. 

For  TL  emission,  where  Eq.  (61)  does  help  determine  E^,,the  effect  of  the 
f^j-function  term  is  again  only  a  small  perturbation  because  x^  is  then 
small.  Thus,  the  f|^  function  is  never  critically  involved  in  the  emission 
current  determination,  anywhere  from  the  SCL  emission  regime  through  to  TL 
emission.  This  verifies  what  was  mentioned  earlier,  that  value  of  f^  is  not 
essential  to  generate  an  SCL-to-TL  transition. 

We  thus  have  shown  that  it  is  the  fy  function  which  critically  determines 
the  emission  current,  while  fj^  primarily  affects  the  final  determination  of 
barrier  position.  We  note  that  this  separation  of  roles  derives  from  the 
algebraic  structure  of  Eqs.  (42a)-(42b).  It  is  this  same  algebraic  structure 
of  these  equations  that  gives  rise  to  the  SCL-to-TL  transition,  rather  than 


the  SCL-to-TL  transition  being  embodied  tflthin  the  explicit  functional  form  of 
the  fy  and  f^^  functions.  The  same  function  for  fy  [Eq.  (55)]  carries  over 
Into  both  SCL  and  TL  regimes.  And  In  both  regimes,  f^^  does  not  significantly 
affect  the  emission  current.  These  are  some  of  the  reasons  why  the  effective 
potential  method,  as  outlined  in  this  report,  can  give  good  approximations  to 
cathode  I-V  and  1-T  curves,  using  just  a  simple  power-law  approximation  for 
the  nonlinear  part  of  the  original  Poisson's  equation  (from  where  fy  and  f|^ 
derive). 

B.  A  MODIFICATION  TO  SCHOTTKY  BARRIER  LOWERING 

The  SCL  emission  current  was  found  to  approximately  be  given  by  "  Jq 

(Child's  Law  Current).  We  next  derive  an  expression  for  In  TL  emission 

limit. 

For  the  Tl.  emission  limit.  It  Is  convenient  to  rewrite  Eq.  (59)  as: 


e^d 


(1  -r^) 


(62) 


Substituting  this  into  Eq.  (61)  and  Ignoring  the  fj^  term  as  being  small  for  TL 
operation,  then  gives: 

1/2 
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(63) 


for  the  Schottky  barrier  lowering  in  TL  emission,  where  Jjiet  ^  "Ic*  Equation 
(63)  agrees  with  the  standard  value  of  the  Schottky  barrier  lowering,  with  a 
new  correction  term  of  order  (Jnet^-lc^*  correction  term  makes  the 

collected  current  less,  due  to  a  persistence  of  some  residual  space-charge 
effects,  at  high  applied  fields. 

The  approximate  solution  to  the  effective  potential  method  equations  thus 
shows  that  ■  J(j  (Child's  Law)  is  valid  over  much  of  the  SCL-reglme,  but  a 

modification  to  Schottky  barrier  lowering  Is  required  for  hlgh-voltages,  due 
to  residual  space-charge  effects  (Eq.  (63)].  The  Child's  Law  limit  Is  valid 
for  barrier  heights  E  >  0,  and  E  vanishes  at  ■  Jj^.  The  two  solutions 


for  low  and  high  fields  thus  Join  together.  In  addition,  the  derivative  also 
can  be  shown  to  be  continuous.  This  property  allows  the  essential  aspects  of 
a  SCL-to-TL  transition  to  be  summarized  in  a  single  equation: 

Jq.  <Jr>  Jg)  (64a) 

Jr,  (Jr-Jc)  (64b) 

•^R  If  <-'r  <  •'c> 

o  c 

to  approximate  the  emission  current  with  both  Schottky  barrier  lowering  and 
space-charge  effects  included.  Equation  (64)  is  Child's  Law  with  a  modified 
Schottky  barrier  lowering. 

In  Fig.  2,  we  examine  how  the  approximate  answer  of  Eq.  (64)  compares 
with  using  the  full  transformation  of  Eqs.  (42a)-(42b).  The  comparison  of 
Fig.  2  was  generated  using  a  work  function  of  2eV,  lOSO’C  temperature,  and  an 
anode-cathode  distance  of  0.0256  cm. 

Figure  2  shows  that  from  SCL-emlsslon  near  the  Richardson  current,  up 
through  the  entire  TL-reglon,  the  approximate  solution,  Eq.  (64),  is 
practically  indistinguishable  from  the  effective  potential  method  calculation 
using  the  full  equations  of  Eqs.  (42a)-(42b).  The  full  equations  do  show  a 
difference  in  the  deep  SCL-reglon  due  to  deviations  from  Child's  Law  in  that 
regime. 

Equation  (64)  thus  represents  a  major  result  of  the  effective  potential 
method.  This  method  has  analytically  calculated  the  approximate  emission 
current  expected  with  both  Schottky-barrier  lowering  and  space-charge  effects 
Included,  summarizing  the  difficult  SCL-to-TL  transition  in  terms  of  a  single 
modified  Schottky-barrier  lowering  term. 

In  the  next  sections,  we  show  how  this  effective  potential  method  com¬ 
pares  to  other  results  for  I-V  curves. 
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Predictions  of  the  "effective  potential"  approxiaation.  A  closed- 
form  equation  which  approximates  the  current  density  vs.  voltage 
[open  circles,  Eq.  (64)  in  text]  is  compared  with  the  more  complete 
transformation  equations  of  the  effective  potential  method  (solid 
line).  The  closed  form  approximation  is  seen  to  be  very  accurate  in 
describing  the  predictions  of  the  present  method  in  the  temperature- 
limited  (XL)  regime.  The  deviations  at  low  voltages  are  due  to 
Child’s  Law  being  Inaccurate  in  the  deep  space-charge  limited  (SCL) 
regime.  A  dashed  line  marks  the  Richardson  current,  which  is  the  SCL 
to  TL  transition  point.  The  calculations  were  done  using  a  work 
function  of  4  >  2.0  eV,  a  temperature  of  1050'’ C,  and  an  anode-cathode 
spacing  of  0.0256  cm  for  all  figures. 
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COMPARISON  TO  SCHOTTKY  AND  CHILD'S  LAW  CURRENTS 


In  this  section,  approximate  cathode  I-V  curves,  generated  by  the  effec¬ 
tive  potential  method,  are  compared  to  Child's  Law  and  the  Schottky  barrier 
lowering  currents.  Here,  the  full  transformations  [Eqs.  (42a)-(42b)]  were 
used  for  and  x^.  Also,  Eq.  (A-4),  in  the  appendix.  Is  used  for  the  f^ 
function,  and  the  large  argument  form  [Eq.  (55)]  is  used  for  fy.  The  fy 
function  critically  determines  the  emission  current,  while  f^  primarily 
affects  the  final  determination  of  barrier  position. 

Figure  3  shows  the  results  of  a  comparison  of  the  effective  potential 
method  to  Child's  Law  and  the  standard  form  of  Schottky-barrier  lo%iering. 
Noteworthy  is  the  fact  that  the  effective-potential  method  predicts  a  net 
current  density  that  is  significantly  below  the  pure  Schottky-llne.  As  noted 
earlier,  this  corresponds  to  residual  space-charge  effects  even  in  the  high 
voltage  limit.  And,  as  was  shown  in  Fig.  2,  Eq.  (64)  is  an  excellent  approxi¬ 
mation  to  the  results  in  this  regime. 

In  Fig.  4,  we  plot  the  same  results  as  Fig.  3,  using  logarithmic  axes. 
This  plot  emphasizes  the  stn'''ture  of  the  I-V  curve  in  deep  SCL-emlssion. 
Child's  Law  is  a  straight  line  of  3/2-8lope  on  this  graph,  and  Fig.  4  shows 
that  significant  deviations  from  Child's  Law  occur  in  the  very  low  voltage 
regime.  This  effect  is  due  to  the  Maxwellian  distribution  of  electrons  being 
explicitly  considered  here,  with  some  high  energy  electrons  always  capable  of 
surmounting  the  energy  barrier. 

In  Fig.  5,  we  plot  the  predicted  barrier  position.  As  noted  earlier,  the 
effective  potential  method  is  expected  to  give  greater  inaccuracies  in  barrier 
position  determination,  while  maintaining  a  small  error  in  net  emission 
current  determination.  The  barrier  position  inaccuracies  should  be  greatest 
in  the  SCL  region,  with  less  of  an  error  in  the  transition  and  TL  regions. 
Figure  5  also  shows  a  dramatic  transition  point,  where  large  changes  in  the 
barrier  position  occur  over  very  small  changes  in  applied  voltage.  This  sharp 
change  reflects  an  inherent  transition  in  the  system  from  SCL-to-TL  emission, 
near  the  Richardson  current  point. 
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Comparison  of  the  present  method  (J[NET])  to  the  Schottky  (J[TLl}  and 
Child  Law  (J[SCL])  limits.  The  net  current  at  low  voltages  is  larger 
than  predicted  on  the  basis  of  Child's  Law  alone.  At  high  voltages, 
the  net  current  also  remains  a  fraction  below  the  Schottky  line,  due 
to  the  persistence  of  space-charge  effects,  deep  into  the 
temperature-limited  regime. 


36 


IJIA/cm^)) 


Fig.  4:  A  plot  of  current  density  vs.  voltage  on  logarithmic  axes.  In  this 
plot,  Child's  Law  is.  a  line  of  slope  3/2.  The  present  method 
predicts  significant  deviations  from  Child's  Law  at  low  voltages  due 
to  the  Maxwellian  distribution  in  energy  of  the  emitted  electrons 
being  taken  into  account  here. 
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Predictions  for  barrier  position  vs.  voltage.  There  Is  a  large 
change  In  barrier  position  within  a  narrow  range  of  voltage.  T 
point  corresponds  to  the  transition  between  SCL  and  TL  emission 
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COMPARISON  WITH  OTHER  APPROXIMATE  METHODS 


In  this  section,  the  approximate  cathode  I-V  curves,  generated  by  the 
effective  potential  method,  are  compared  to  two  prior  methods,  one  by  Scott, 
and  which  Rasker^^  has  analytically  approximated,  one  by  Longo.^^ 

Historically,  Langmuir^  modified  the  simple  power  law  of  Child^  for 

emission  current  density*  As  a  function  of  applied  voltage,  the  Langmuir 
solution  begins  in  the  almost-retardlng-field  (ARP)  regibe  with  significant 
deviations  from  Child's  law,  follomd  by  a  range  where  Child's  law  is 
approximately  valid,  then  it  abruptly  bends  over  to  the  Richardson  current 
value,  first  reaching  the  Richardson  current  value  at  some  critical  voltage, 

Vl* 

Schottky^  published  his  modification  to  the  Richardson  current  the  same 
year,  but  adding  Schottky  barrier  lowering  to  the  Langmuir  solution  has  been 
slow  to  develop.  Scott^^  attacked  the  problem  by  assuming  the  Langmuir 
solution  is  entirely  valid  for  voltages  less  than  V^,  even  in  the  presence  of 
Schottky  barrier  lowering.  For  voltages  greater  than  he  replaces  the 
applied  field  in  the  Schottky  barrier  lowering  effect,  with  an  "effective 
field"  derived  from  the  local  Langmuir  potential.  Using  this  technique,  he 
circumvents  the  inhomogeneous  part  of  the  nonlinear  Poisson's  equation*  The 
joining  of  solutions  at  was  noted  by  Scott  to  give  rise  to  a  "kink"  in  the 
solution,  where  the  current  vs.  voltage  goes  from  horizontal  (zero  slope)  to  a 
substantial  slope  just  after  V^. 

Longo's  method is  almost  seml~empirical.  He  finds  that  some  data  on 
actual  cathodes  seem  to  fit  the  equation: 

(1/Jnet>  "  +  (1/J[TL]),  (65) 

where  J[TL]  is  the  standard  temperature-limited  value  for  current  density  with 
Schottky  barrier  lowering  Included,  and  where  is  Child's  law.  The  appeal 
of  Longo's  equation  is  its  simplicity* 

In  some  sense,  the  present  work  is  similar  to  both  Scott  and  Longo.  Like 
Scott,  we  circumvent  dealing  with  the  full  inhomogeneous  nonlinear  Poisson's 
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equation,  making  both  Scott  and  this  solution  different  approximate  answers  to 
the  same  problem.  Scott  uses  the  Langmuir  equation  for  the  nonlinearity  part, 
"patchlng^in"  Schottky  barrier  lowering  using  an  effective  field.  Here,  we 
dealt  explicitly  with  the  inhomogeneous  Poisson's  equation  which  describes 
Schottky  barrier  lowering,  and  we  "patched-in”  the  nonlinearity  by  using  an 
effective  potential. 

However,  like  Longo,  we  sought  to  preserve  as  much  simplicity  as  possible 
In  our  approximate  solution,  by  keeping  our  solution  close  to  the  domain  of 
elementary  functions.  The  hope.  Implicit  In  Longo 's  work.  Is  to  have  a  single 
equation  describing  cathode  behavior  from  the  space-charge  limited  (SCL)  to 
the  temperature-limited  (TL)  regime.  That  goal  Is  realised,  to  a  much  greater 
degree  of  accuracy  than  Longo 's  equation,  by  Bq.  (64)  of  this  work. 

Figure  6  plots  the  current  vs.  voltage,  using  logarithmic  axes,  for  the 
same  conditions  as  In  Fig.  4:  1050*C,  2  eV  work  function,  and  an  anode-cathode 
distance  of  0.02S6  cm.  It  compares  the  result  of  this  work  (solid  line)  to 
the  prior  work  of  both  Scott  and  Longo  (dotted  lines). 

At  low  voltages,  Longo's  equation  asymptotically  approaches  Child's  Law 
from  below,  with  the  other  methods  giving  more  current  than  Child's  Law  due  to 
contributions  from  the  high  energy  tall  of  thermally  excited  electrons  in  the 
emitter.  The  present  work  overestimates  the  current  In  this  low-voltage  ARF 
regime  the  most,  mainly  due  to  neglecting  the  explicit  contribution  of 
electrons  which  can  form  a  trapped-cloud  near  the  emitter  surface.  Scott's 
work  uses  Langmuir's  solution  at  these  lowest  voltages,  but  even  at  eV  values 
comparable  to  kT,  this  work  Is  accurate  to  within  a  factor  of  two. 

On  the  broad  scale  of  Fig.  6,  covering  more  than  four  decades  in  both 
voltage  and  current,  Scott  and  this  solution  are  similar.  Longo's  equation  Is 
low  by  a  fair  margin,  except  In  the  deep  temperature-limited  region. 

The  critical  transition  or  knee  region  between  SCL  and  TL  operation  Is 
shown  on  an  expanded  scale  In  Fig.  7.  In  this  transition  region,  the  general 
shape  predicted  by  both  Scott  and  this  work  Is  similar.  Scott  gives  a  net 
emission  current  slightly  higher  by  a  few  percent  than  this  work,  although  an 
exact  numerical  calculation  Is  needed  to  determine  which  method  Is  relatively 
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Fig.  6:  Comparison  of  current  density  vs.  voltage  predictions  among  several 
methods.  The  approximate  methods  of  this  work,  Scott,  and  Longo  are 
compared  over  several  orders  of  magnitude  in  both  current  and 
voltage.  The  Longo  equation  is  low  by  a  fair  margain  except  in  the 
deep  temperature-limited  region,  with  the  other  methods  giving 
similar  results  on  this  scale. 
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Comparison  as  in  Fig.  6,  on  an  expanded  scale,  emphasizing  the  “knee*' 
or  transition  .region  from  SCL  to  TL  operation.  Child’s  Law  and  the 
Schottky  barrier  lowering  limit  are  also  shown  for  comparison.  This 
work  predicts  values  that  are  a  few  percent  lower  In  the  knee  region 
than  the  Scott  approximation. 


more  accurate  here.  Both  break  away  from  Child's  Law,  giving  more  SCL  cur¬ 
rent,  with  this  work  being  more  gradual  than  Scott.  In  the  TL  region,  the 
Scott  result  seems  to  approach  the  Schottky  line  faster  than  this  work. 

In  Fig.  8,  we  plot  the  relative  error  among  the  several  methods.  This 
work  was  chosen  as  the  0%  baseline.  The  vertical  axis  is  percent  error  in 
emission  current,  spanning  only  2SZ  net  difference,  while  the  horizontal  axis 
spans  four  orders  of  magnitude  in  voltage. 

Figure  8  clearly  shows  that  this  work  remains  bounded  away  from  the  pure 
temperature-limited  line  (J[TLl)  of  Schottky  barrier  lowering.  [Equation  (64] 
is  also  plotted  on  Fig.  8.  It  is  a  one-line  approximation  to  this  work,  and 
it  has  the  accuracy  of  Child's  law  in  the  SCL  regime. 

Examining  the  TL-region  of  the  Scott  solution  and  the  Longo  equation  in 
Fig.  8  shorn  that  they  both  approach  the  J[TL]  line  too  quickly.  This  effect 
is  due  to  both  methods  not  taking  into  account  the  residual  space-charge 
effects  in  the  TL-reglon. 

The  comparison  of  methods  in  Figs.  6-8  supports  the  usefulness  and  over¬ 
all  accuracy  of  this  new  method.  In  addition,  the  closed-form  analytic  cor¬ 
rection  to  Schottky  barrier  lowering  Eq.  (64)  retains  practically  all  of  the 
essential  features  associated  with  knee-roundedness,  and  it  can  be  used  as  a 
simple  approximation  to  the  TL-emlsslon  region  for  uniform  work  function 
emitters. 


Fig.  8:  Comparison  of  relative  percent-error  among  the  different  methods. 

Both  this  work  and  Scott  are  similar  although  this  work  tends  to  be 
smoother.  The  closed-form  approximation  of  Eq.  (64).  which  Is  a  new 
correction  to  Schottky  barrier  lowering,  Is  also  shown.  Equation 
(64)  has  the  accuracy  of  Child *s  Law  In  the  SCL  regime,  and  virtually 
the  same  accuracy  as  this  work  In  the  TL  region.  The  high  voltage 
region  also  shows  that  the  Scott  and  Longo  methods  approach  the  usual 
Schottky  barrier  lowering  current  J[TL]  too  quickly. 


V.  CONCLUSIONS 


In  this  work,  ve  have  developed  a  new  technique  for  approximating  the 
emission  current  from  a  planar  thermionic  emitter  In  a  diode-type  geometry. 

In  general,  this  procedure  reduced  the  nonlinear  Poisson's  equation  to  finding 
the  roots  of  two  algebraic  equations.  These  two  roots,  labelled  and  x^, 

correspond  to  the  net  barrier  height  for  electron  emission,  and  net  barrier 
distance  from  the  cathode. 

This  method  was  then  applied  to  the  problem  of  a  uniform  work  function 
surface  with  both  Schottky  barrier  lowering  and  space-charge  effects 
Included.  The  equations  for  the  roots  and  x^  In  this  case  were  solved 
using  root-finding  routines,  and  the  net  emission  current  was  compared  to 
previously  published  approximations^^* for  the  uniform  work  function  case. 
Scott's  numerical  method  and  this  analytic  one  were  found  to  quantitatively 
agree  on  cathode  emission  current  to  within  •>  lOZ  over  a  wide  range  of  applied 
voltage. 

In  addition,  the  analytic  formalism,  although  approximate,  allowed  us  to 
derive  a  new  closed  form,  representing  a  first  correction  to  Schottky  barrier 
lowering,  due  to  space-charge  effects  being  Included.  Physically,  this  new 
correction  term,  not  Included  In  any  previous  approximation,  corresponds  to  a 
persistence  of  space-charge  effects  deep  Into  the  temperature-limited 
regime.  This  closed  form  specifies  the  roundedness  of  the  knee-region  between 
space-charge  limited  (SCL)  and  temperature-limited  (TL)  emission  for  the 
uniform  work  function  case.  It  determines  the  SCL-to-TL  transition  point  as 
being  the  Richardson's  current.  The  closed  form,  approximating  the  emission 
current  density  for  a  uniform  work  function  surface,  was  found  to  be: 

Jnet  “  ^  (66a) 

■'n.t  •  -’ll  *  ‘"o  “  -  ■’C  >  Jr- 


(66b) 


In  this  equation,  •^net  is  the  total  eiBlsslon  current  density,  J|^  and  are 

the  Richardson  and  Child  Law  current  densities,  and  M  Is  the  usual  value  for 

o 

Schottky  barrier  lowering.  This  simple  approximation  for  Is  also  contin¬ 

uous  at  Jj^,  with  a  continuous  first  derlvltlve. 
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APPENDIX;  THE  fy  AND  FUNCTIONS 


In  this  appendix,  we  solve  the  two  integrals,  Eqs.  (39a)  and  (39b), 
determine  the  nonlinear  fy  and  terms. 

A.  THE  FUNCTION 

The  Integral  representing  the  term  is: 


ft (A)  -  (1  +  u)  du 


(4/0 


1/2 


u— 1 


Defining: 


-  It  A/4, 


The  fjj-f unction  can  then  be  written  as: 

R 

f,^(A)  -  j  (1  -  I  )  dr/I(r^  +  1)^^^  +r] . 
These  Integrals  can  be  done  exactly,  resulting  in: 

-  -Jr  An  {R  +  (R^  +1)^^^]  - 

(3R^’-4R) _ 

6(R^-2)  +  6(R^-2)  (R^  +1)^^^ 


for  the  ft  function. 

.  THE  fy-FUNCTION 

The  integral  representing  the  fy-term  is: 


which 


(A-1) 


(A-2) 


(A-3) 


(A-4) 


6 


Using  the  substitution: 


r 


2 


(.A/4)  s 


the  fy  function  can  then  be  written  as: 


(A-6) 


f  (A)  -  [l-({  dr/[(r2  +  l)  ^  ^ 

y  r-0  * 


The  term: 


E(R)  =  / 


R 


r-0  r  +  [r^ 


(A-8) 


Is  an  elliptical  Integral,  the  other  term  being  an  elementary  function. 

Since  R^  «  (eV^/kT),  where  la  the  applied  accelerating  voltage,  one  Is 
generally  Interested  In  the  large  R  limit  for  the  f^  function.  A  large  R 
expansion  for  the  elliptical  Integral,  E(R),  Is: 


E(R)  «  R^^^-  C  +  (12R^^^)"^  +  {Order 
o 


(A-9) 


where  Is  an  Integration  constant. 

The  particular  value  of  C^,  specified  by  choosing  the  lower  limit  of  Eq. 
(A-8)  to  be  r-0,  is: 


C  «  0.677771,  (A-10) 
o 

determined  numerically.  It  Is  Important  to  note  that  Is  a  parameter  In  the 
large  R  expansion  of  E(R)  which  Is  sensitive  to  the  behavior  at  small  R.  Thus 
Cg  Is  sensitive  to  details  of  the  effective  potential  which  was  chosen. 


The  ocher  Integral  appearing  in  the  fy  function: 


/  r  dr 

F(R)  “  j  2  1/2 

r-0  r  +  (r^  + 

can  be  done  exactly  using  Che  substitution  r  -  sinh(z). 
for  F(R)  is: 


(A-11) 


The  large  R  expansion 


P(R)  -  R^/4  -  Iln(R)  +  Cj  1/8  -  (32  R^)"^  +  [Order  (R“^)].  (A-12) 

is  another  constant  of  integration  in  the  large  R  expansion  which  is  sensi¬ 
tive  to  the  details  of  the  effective  potential »  similar  to  The  particular 
value  of  Cj  for  Eq.  (A-10)  is: 


Cj  -  (ln2  -  1/4)  «  0.443147. 


(A-13) 


This  gives  the  following  large  R  expansion  for  the  function  fy: 


■7T37r  * 


t31n(R)  +  3Cj  +  2] 


-  [Order  (R"^)], 


(A-14) 


8R  2  R 


where : 


R^  -  (itA/4) 


(A-15) 


The  important  conclusion  for  Che  fy  function  is  that  in  the  region  of  interest 
(large  R),  fy  has  an  expansion  whose  leading  term  is  (9/8R).  Also,  via  Cq.  it 
has  a  first  correction  term  whose  magnitude  depends  on  Che  details  of  Che 
potential  between  the  electrodes,  from  the  collector  all  the  way  back  to  Che 
near  cathode  region. 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  Is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  r^w  military  space  systems.  Versatility 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  person¬ 
nel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  systems.  Expertise  in  the  latest  scientific  developments  is 
vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The  labora¬ 
tories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  environmental  hazards,  trace  detection;  spacecraft  structural 
mechanics,  contamination,  thermal  and  structural  control;  high  temperature 
thermomechanics,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  development 
Including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control, 
atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency 
standards,  and  environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence  and 
microelectronics  applications. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  lidar,  and  electro¬ 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  imaging;  millimeter  wave,  microwave  technology, 
and  RF  systems  research. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospberlc,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy, 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 
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